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We study theoretically the effects of spatial confinement on the phoretic motion of a dissolved 
particle driven by composition gradients generated by chemical reactions of its solvent, which are 
active only on certain parts of the particle surface. We show that the presence of confining walls 
increases in a similar way both the composition gradients and the viscous friction, and the overall 
result of these competing effects is an increase in the phoretic velocity of the particle. For the case 
of steric repulsion only between the particle and the product molecules of the chemical reactions, 
the absolute value of the velocity remains nonetheless rather small. 



PACS numbers: 89.20.-a, 07.10.Cm, 82.56.Lz 

1. INTRODUCTION 

Recent years have witnessed a growing technological, 
experimental, and theoretical interest in scaling standard 
machinery down to micro- and nano-scales as needed for 
the development of "lab on a chip" devices. For appli- 
cations in, e.g., drug-delivery systems or micromcchanics 
one of the most challenging problems at this stage is to 
develop ways to enable small-scale objects to perform au- 
tonomous, controlled motion^i^. Although the research 
in this area is still in its early stages, several such pro- 
posals have already been tested experimentally (see, e.g., 
Refs. [313110) • A review of the recent progress in this 
field can be found in Ref. 0. 

Whitesidcs and co-workers proposed a design of self- 
propelling devices based on an asymmetric decoration 
of the surface of small objects by catalytic, active sites 
promoting a chemical reaction in the surrounding liq- 
uid medium^. This asymmetric distribution can pro- 
vide motility through a variety of mechanisms, such as 
surface tension gradients and/or cyclic adsorption and 
desorptioi^ii. The use of an asymmetric surface distri- 
bution of a catalyst has been further proposed for an au- 
tonomous diffusiophoretic motion emerging as a result of 
self-created density gradients^. An experimental realiza- 
tion of such systems, using platinum coated polystyrene 
spheres, has been recently reported^. The issue of de- 
signing optimal distributions of a catalyst for spherical 
and rod-like particles has also been approached^. 

For most of the applications in biological systems or in 
'lab on a chip'-type devices one has to deal with a compli- 
cated internal structure of the systems, such as networks 
of narrow channels or pores and various impenetrable 
impurities. In some situations, one may even encounter 



a quasi two- or one-dimensional behavior, like in the 
cases of bacteria motion on planar nutrient substrates^ 
and of the motion of small particles within a biologi- 
cal membranei^ or in a monolayer adsorbed on a three- 
dimensional {3d) liquid subphase>ii. Thus spatial con- 
finement is a relevant feature so that the assumption of 
the presence of an unconfined 3d bulk reactive solveni^i*^!^ 
may break down. Intuitively, one expects that spatial 
confinement docs influence the resulting motion of self- 
propelling objects like the ones discussed above. But a 
priori it is not clear if such effects are significant. 

Starting from the model used in Ref. here we study 
the effects of spatial confinement on the phoretic motion 
of a particle that generates number density gradients of 
the product molecules emerging from the chemical reac- 
tions in the depleting Newtonian liquid solvent. We shall 
focus on the simple case in which the particle and the 3d 
solution of solvent and product molecules are bounded by 
a spherical shell because (i) this simple geometry allows 
for an exact solution and (ii) it has been experimentally 
shown for a variety of phoretic systems that confinement 
effects are dominated by the smallest confining length 
scalei^. Thus this geometry has a paradigmatic charac- 
ter. 

Similarly to the earlier studiesi'^'^'^'^, our work is 
based on adopting the standard theory of phoresis for 
the present case, in which the gradients are self-generated 
rather than being produced and maintained by external 
sources. In doing so, one is bound to make a number of 
assumptions that are either already present in the clas- 
sical theory of phoresis, or arise as a result of mapping 
the description of such "active" surface particles onto the 
framework of a theory developed to describe the case 
of inert particles immersed in a pre-defined, externally 
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controlled concentration gradient. Since in the literature 
these assumptions are often overlooked or not spelled out 
explicitly, we consider it as necessary to provide also a 
critical discussion of the significant assumptions involved 
by this mapping, as well as of some of those assumptions 
implicitly contained in the standard theory of phoresis. 
Accordingly, the outline of this paper is as follows. In 
Section [2] we define the model and discuss some general 
aspects of systems with self-generated motion, with par- 
ticular emphasis on the assumptions involved in adopting 
the standard theory of phoresis. Section [3] is devoted to 
the computation of the diffusiophoretic velocity. This 
includes the calculation of the diffusiophoretic slip ve- 
locity and the determination of the phoretic hydrody- 
namic flows and density profiles of the product molecules 
around a self-propelling particle. The results are dis- 
cussed in Section |4l and we conclude in Section [5] with 
a brief summary of our results and general conclusions. 
Important details of our calculations are presented in the 
Appendices El IHl andO 



2. THE MODEL 

2.1. General aspects. 

The system we consider is shown in Fig. [1] It con- 
sists of an impermeable, spherical particle of radius R 
with a point-like catalytic site (black dot in Fig. [T|) on 
its surface, which promotes the chemical conversion of 
a surrounding solvent into product molecules of diame- 
ter a (small hatched circles in Fig. [T|). The particle and 
the surrounding solution (solvent plus the solute, i.e., the 
product molecules) of viscosity fj, are enclosed in a con- 
centric, impermeable, spherical shell of radius Ri ~ -qR 
iv > !)■ 




FIG. 1: An impermeable, spherical particle of radius R 
with a point-like catalytic site at r^ = —Re^ (depicted as 
a black dot) on its surface, enclosed by a concentric, imper- 
meable, spherical wall of radius _Ri = rjR {rj > 1). The prod- 
uct molecules of diameter a are shown as small horizontally 
hatched circles. 



In general, such a chemical conversion of the solvent 
gives rise to several types of product molecules. Here 
we shall focus on the particular case in which the chem- 
ical conversion of a solvent molecule (A) leads to two 
molecules only (A' and B), one very similar in size and 
properties with the solvent itself (A' w ^4), the other one 

(B) different [A ™* A' + B]; in the following this latter 
is denoted as "product molecule" and plays the role of 
a solute in the solvent. In other words, we consider a 
situation in which the net result of the chemical conver- 
sion can be approximated as the generation of a solute, 
solely, and the reaction does not lead to a solvent deple- 
tion near the catalytic site which otherwise would reduce 
the current of solvent towards the catalytic site acting as 
a solvent sink. For example, this is approximately the 
case for the Ft catalyzed decomposition of hydrogen per- 
oxide (II2O2) in aqueous solution into water (H2O) and 
oxygen (O2) molecules, as discussed in Refs. [1100. In 
these experimental studies the oxygen plays the role of 
the product molecule the properties of which differ sig- 
nificantly from those of the solvent. Here the solvent is 
actually a binary liquid mixture of H2O and H2O2 for 
which H2O is passive and does not participate in the 
chemical conversion. We note that actually it is rather 
difficult to assess whether or not H2O and H2O2 can be 
treated as being the same in such systems, and proba- 
bly the answer would be on a case by case basis. This 
is not only because the H20-particle and H202-particle 
interactions may be quite different for different materials 
(the simplest example is exactly the catalytic decompo- 
sition of H2O2 on Ft), but also because of the sensitivity 
of phoresis to the details of the solvent mediated interac- 
tions - and the weakly acidic nature of H2O2 may play 
an important role here. In providing this example as a 

possible realization of a catalyzed reaction A A + B 
we had to rely on: (i) the statement in Ref. [second en- 
try, bottom of left column on page 13427 therein], which 
is based on some previously published results by Fhibbs 
and GigucroH^ that hydrogen peroxide and water in con- 
tact with Au have almost identical interfacial tensions 
(this translates into similar particle-solvent interactions), 
and thus that there is no additional solvent-density gra- 
dient to be considered, and (ii) on the implicit statement 
in the same Ref. [second entry, Eq. (8) therein] that 
the effective interaction of the O2 molecule with the Au 
surface of the rod does not depend on the composition of 
the solvent H2O - H2O2 mixture. 

We thus assume that the reaction at the catalytic site, 
(located at Vg = —Rsz, where is the unit vector of the 
2;-axis) , acts effectively only as a point-like source of prod- 
uct molecules^ of diameter a which are diffusing in the 
solvent with diffusion coefficient D. In passing we note 
that this is in contrast to the situation considered in Ref. 
0, where the catalytic site acts both as a source for the 
product and as a sink for the solvent, i.e., the generation 
of a product molecule is accompanied by the annihilation 
of a solvent one so that A"^ B which implies that B can 
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only be a different configuration of because it must 
have the same chemical constituents. In the situation 
considered here, the amplitude of the production rate of 
B is denoted as B{t). We shall neglect any interaction be- 
tween the product molecules. Thus the number density of 
product molecules is considered to be so low that among 
themselves they behave like an ideal gas. There is an 
interaction potential between the product molecules and 
the moving particle (see Appcndix[X|, which includes the 
impermeability condition at the surface of particle. The 
interactions between the product molecules and the sol- 
vent are accounted for in an effective way via the Stokes 
- Einstein expression D = ksT / {i^n^a) for the diffusion 
coefficient D of the product moleculesi^, where fcs is the 
Boltzmann constant and T is the temperature. 



2.2. Discussion of the model in the context of the 
standard theory of phoresis. 

The presence of a source of solute (product molecules) 
on the surface of the particle creates a non-uniform and 
time dependent distribution of solute in the solution (see 
Fig. [ij, i.e., a non- uniform composition of the solution. 
Because the understanding of the way in which such a 
non-uniform composition gives rise to phoretic motion is 
not straightforward, we discuss here in some detail how 
the general model introduced in Section 12.11 can be put 
into the context of the standard theory of phoresis^. A 
recent clear exposition of the basic concepts and a dis- 
cussion of different phoretic transport scenarios can be 
also found in Ref. [H. 

If the production rate of the source is not very large 
and the diffusion coefficient of the product molecules 
through the solvent is not too small, the solute density 
varies smoothly in space and slowly in time. This justi- 
fies the assumption of local equilibrium and the defini- 
tion of a position- and time-dependent chemical potential 
of the solution (the spatial gradients of which describe 
the diffusion of the solute away from the source). In a 
quasi steady-state of the solution, corresponding to slow 
time- variations of the composition, the net solute current 
is approximately zero and the solute density gradients 
are balanced by pressure gradients [see also Eq. (B3) in 
Ref. [H. Because in the absence of external body forces 
the mechanical equilibrium for a liquid solution is gener- 
ally established much faster than the chemical one (see, 
e.g., Ref. HI), it is physically plausible to assume that 
far from the boundaries (i.e., far from the wall and far 
from the particle surface) the pressure adjusts instan- 
taneously to accommodate the spatially varying solute 
density profilo^^ and keeps varying slowly with time (on 
time scales much larger than the diffusion time ~ / D) 
because the overall density of the solution increases in 
time. In other words, the pressure is determined by the 
solute density profile and is obtained from the equation of 
state of the solvent-solute binary mixture considered as 
a thermodynamic system in local equilibrium described 



by the corresponding free energy density. Clearly, the 
situation is different in the region of contact between 
the solution and the particle surface (and also of con- 
tact with the wall) where the interactions of the particle 
with the product molecules (solute) and with the solvent 
molecules are relevant. This interaction disturbs the dis- 
tribution of solute molecules, net currents of solute and 
solvent result, and the pressure field becomes a quantity 
determined by the incompressibility requirement for the 
hydrodynamic flow, and not by an equation of state. We 
shall discuss this important point below. 

At a given time the spatial variations of the number 
density p of the solute and of the number density Psoiv of 
the solvent are characterized as follows: Psoiv is constant 
throughout the solution apart from the close vicinity of 
the particle surface and of the surface of the confining 
sphere. There psoiv vanishes and reaches its constant 
bulk value via density oscillations which arc induced by 
local packing effects caused by the finite diameter agoiv 
of the solvent particles. Away from any kind of surface 
phase transition the approach of the bulk value occurs 
exponentially on the scale of the bulk correlation length 
isoiv of the solvent. (The presence of long-ranged dis- 
persion forces causes power-law decays; however, their 
amplitude is small and is neglected in the present con- 
text.) Away from bulk phase transitions S,soiv is compa- 
rable with the range of the interaction potential between 
the solvent molecules, which in turn is proportional to 
Qsoiv with a prefactor of order unity. Within this picture 
Psoiv does not vary along the surface of the particle. 

The solute number density p is characterized by two 
important features. On the scale of the system size 
R ^ a, p varies due to the diffusion process. On the 
much smaller length scale of the solute diameter a, p 
also varies near the particle surface and near the surface 
of the confinement as discussed above for the solvent. 
Since, however, according to our earlier assumption the 
solute particles can be considered to form an ideal gas, 
near the walls p varies proportional to exp(— /?'!'), where 

is the effective substrate potential (see Appendix |^ 
in the sense that it describes the interaction of the so- 
lute molecule with the substrate in the presence of the 
solvent, and (3 = l/ZcsT. Typically the range (5 of 
is proportional to the solute diameter a. Accordingly 
the solute molecules interact directly with the particle 
only if they are within a thin surface film of thickness S, 
which is assumed to not be deformed by the motion of 
the particlei^. 

The standard theory of phoresis assumes a, S ^ o-soiv 
so that within this surface film of thickness 5 the solvent 
can be considered as a continuum with constant Psoiv 
We note that this assumption is rather convenient from 
a computational point of view. However, this assumption 
is at odds with the basic underlying model according to 
which the solute particle is created from a solvent parti- 
cle via a catalytic reaction. Under these circumstances 
one expects a w agoiv (see the above example of the reac- 
tion 2 H2O2 ^ 2 H2O + O2). In the absence of a more 
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detailed theory of phoresis, which treats the sizes of the 
solvent and solute particles on the same footing, we pro- 
ceed along the lines of the standard theory. Nonetheless 
we point out that there is an urgent need for improve- 
ment. It might be that due to the fact that continuum 
hydrodynamics remains quantitatively reliable down to 
surprisingly small length scales (here asoiv), it is con- 
ceivable that the aforementioned continuum description 
yields reliable results, too. But this remains unprovcn. 

According to the standard theory a hydrodynamic de- 
scription applies within the aforementioned surface film. 
Within this picture the solute molecules and their effec- 
tive interaction with the particle are replaced by a cor- 
responding distribution of "point forces" acting on the 
solvent in the film. Within the limitations of such an 
approach the hydrodynamic description of the solution 
naturally splits into that of an "inner" region formed by 
the surface film and that of the "outer" region formed 
by the exterior space beyond the surface film. Following 
Rcfs. 8. 12.,17.. the ensuing asymmetric, non-uniform so- 
lute number density p(r, t) around the particle will give 
rise, within the surface film only, to a pressure gradient 
along the surface of the particle. This is because the 
surface film is very thin on the scale of the system size 
R so that the equilibration of the composition profile of 
the solution within the surface film can be assumed to 
be fast along the direction normal to the surface of the 
sclf-propcUing particle compared with the diffusional re- 
laxation time of the composition gradient along its sur- 
face, which typically involves a length scale of the order 
of the particle size R. Therefore, within the surface film 
the solute density in the direction normal to the surface 
of the particle is given by a Boltzmann distribution cor- 
responding to the local equilibrium configuration in the 
presence of the effective interaction potential between the 
particle and the solute molecules, with a prefactor which 
depends on the position along the surface of the particle 
(for details see Appendix [XI which builds on Ref. [T^ . 
Mechanical equilibrium of the solvent within the surface 
film along the direction normal to the surface of the par- 
ticle (no flow along this direction) requires the pressure 
gradient along the normal to be equal to the body force 
densities due to the effective particle-solute interactions. 
Therefore the pressure within the surface film differs from 
the "outer" pressure field by an "osmotic pressure" term, 
i.e., a term proportional to the extra solute density p in 
excess to the constant solvent density Psoiv Since this 
osmotic pressure varies along the surface of the particle, 
there is a gradient of pressure along the surface of the 
particle. This lateral pressure gradient is not balanced 
by any body force and thus generates shear stress within 
the surface film. Therefore it induces hydrodynamic flow 
of the solution along the surface of the particle and en- 
tails motion of the particle with a velocity V(i). Because 
the system has azimuthal symmetry, the motion is along 
the z-axis, i.e., V = Vsz- The hydrodynamic flow of the 
solution, the diffusive transport of the solute, and the 
coupling between the two giving rise to phoretic motion 



of the particle are analyzed in detail in the next section 
(see also the Appendices MUj) . 



3. COMPUTATION OF THE 
DIFFUSIO-PHORETIC VELOCITY 

Our calculation of the diffusiophoretic velocity pro- 
ceeds along the lines of Ref. [13 and is based on dividing 
the problem into an inner one - within the film around 
the particle surface discussed in Subsec. 12.21 and an outer 
one - beyond the range of the effective interaction be- 
tween the solutes and the particle. In view of the argu- 
ments presented in Sec. [21 we base our analysis on the 
following assumptions: 

(i) The chemical reaction leads to a change in the solute 
density only. The thickness of the surface film, defined by 
the range of the effective interaction between the prod- 
uct molecules (solute) and the particle, is much smaller 
than the particle radius R. The spatial variations of the 
solute number density along the surface of the particle 
occur over length scales of the order of R, which allows 
one to use the approximation of a locally planar inter- 
face. 

(ii) The number density of product molecules is suffi- 
ciently low, so that the solute can be viewed as an ideal 
gas and the solvent-solute mixture behaves like an ideal 
dilute solution. In this case, the pressure gradient is sim- 
ply proportional to the gradient of the number density of 
the product particles (see Appendix [K^ . 

(iii) Temporal variations of the number density of the 
product molecules due to their creation occur on time 
scales much longer than those needed for the fluid flow 
(as seen from the moving particle) to relax to a steady 
state corresponding to the number density profile at that 
moment. 

(iv) The flow field of the solution within the surface film 
can be described by the laws of hydrodynamics. 

Additionally, we assume that both the Reynolds num- 
ber Re ~ psoivV R/ p., where psoiv is the mass density 
of the solvent, and the Peclet number Pe ~ VR/D are 
small, such that one can approximate the hydrodynamic 
description with the creeping flow (Stokes) equations and 
disregard the convection of the solute compared to its 
diffusive transport. Here we have assumed that the mag- 
nitude of the hydrodynamic flow u is similar to that of 
the phoretic velocity V; this assumption is supported a 
posteriori by the fact that for our system the phoretic 
velocity is basically the average of the slip-velocity over 
the surface of the particle [see, c.f., Eqs. ([T]) and PU)) ]. 
For a yum size particle moving through water (density 
Psoiv = 10"^ kg/m , viscosity p. = 10~^ Pa s) with a ve- 
locity of the order of pm/s, which is typical for phoresis, 
one has Re ~ 10^^. For the diffusion at room temper- 
ature {kBTroom 4 X IQ-^i J) of O2 (a ~ 10-1° m) in 
H2O2 {p — Pa s), the Stokes-Einstein relation leads 
to an estimate Z? ~ 4 x lQ~^w?ls for the diffusion coef- 
ficient (in agreement with Ref. [j), and thus Pe ~ 10"^. 
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Therefore the latter assumptions arc justified. Note that 
if one uses i?i rather than i? as a characteristic length 
scale, the above results imply that the Re and Pe num- 
bers remain both very small as long as 77 < 10. 



3.1. Diffusiophoretic slip-velocity 

As discussed in Sec. [H the pressure gradient along the 
surface of the particle, induced by its interaction with 
the non-uniform distribution p{r, t) of product molecules, 
leads to flow of the solution relative to the particle. As 
shown in Appendix the hydrodynamic flow within the 
surface fllm translates into a (phoretic) slip-velocity. 



v,(r,t) = -6V,p(r,i), for |r| = i?4 



(1) 



at the outer edge R+ > Roi the surface fllm as a bound- 
ary condition for the hydrodynamic flow in the outer re- 
gion. In this equation denotes the projection of the 
gradient operator onto the tangential planes of the sur- 
face of the particle. 
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where and arc the polar and azimuthal unit vectors, 
respectively, while 



-A 



(3) 



is an effective "mobility", and A = y |A| is a character- 
istic length scale. The latter is determined by the effec- 
tive interaction potential ^ between the particle and the 
product moleculesi^ [which determines their distribution 
within the surface fllm along the direction y normal to 
the particle surface within the local coordinate system 
(Fig. [3]) or the radial direction here]: 



A = / dyy 



(4) 



In the case of a purely steric repulsive interaction, one 
hasi^ A = -a^/8 so that A = 0/(2^2). For a single 
catalytically active site as shown in Fig. [1] p and thus v,, 
do not depend on the azimuthal angle (p. 



3.2. Phoretic hydrodynamic flow 

We focus here on the quasi-static approximation (iii) 
described above and the limit of low Re numbers. In 
the laboratory frame, in which the particle moves with 
a yet unknown velocity 'V{t]irj), the hydrodynamic flow 
u = u(r; t, rf) at time t in the domain beyond the surface 
fllm around the particle obeys the steady-state force-free 



Stokes equations 

V • n = 0, V • u = 0, i?+ < |r| < i?i , (5) 

where II = —pi + /iS is the pressure tensor, p is the hy- 
drostatic pressure, and S is the shear stress tensor, i.e., 
Sc(/3 = dua/dxp -\- dup/dxa, subject to boundary con- 
ditions (BC) of no-slip at the wall at Ri and prescribed 
slip velocity on the surface i.e.. 



V(t;,7) 



|r| = i?.i 



0. 



(6) 



Note that the parametric dependences of the flow fleld 
u on t and rj stem from the boundary conditions. Note 
that because for the outer problem the variations of the 
flow fleld and those of the number density of the prod- 
uct molecules are over length scales that are much larger 
than (5, one can replace everywhere in the calculations 
the (unknown) value of i?+ by i?; wc shall use this ap- 
proximation in the following. 

The solution of the Stokes equations with boundary 
conditions on spherical surfaces is based on expressing 
both the flow fleld u and the pressure fleld p as series 
of solid harmonic&i^ii^ Ki (r, 9) {I S Z) , which are the 
eigenfunctions of the 3d Laplace operator. In Appendix 
IB] we provide a brief outline of the general method for 
solving Eq. ([5|) in spherical coordinates and derive the 
solution obeying the BCs given by Eq. [HI 

Equation ([6]) at Ri corresponds to a no-slip condition 
imposed on a wall flxed with respect to the laboratory 
frame. This deserves further discussion. In Eq. ([6]) the 
BC is that u at the surface of the moving particle, i.e., 
u||r-Rp(t)|=fi+ where Rp(t) = J^Y(t')dt' is the position 
of the particle center, takes the value V(i; ri)+Vs. As long 
as |Rp(<)| ^ concentricity holds and one obtains the 
flrst of the two BCs given by Eq. If, however, the 

velocity V is not small, we readjust the center of Ri in 
order to impose concentricity using the following "pro- 
tocol" : the particle is allowed to move for a short time 
St with the instantaneous velocity V(t) while the spher- 
ical shell at Ri is flxed, after which the spherical shell 
is displaced by 'V{t)6t to its new flxed position in such 
a way that it does not perturb signiflcantly the density 
and flow flclds; this procedure is then repeated. It is not 
clear to which extent this latter assumption, which al- 
lows us to obtain an analytical solution, can be realized 
experimentally. But it is expected that the results we 
derive for the present geometrical setup are relevant also 
for more general geometries^S^ such as a particle moving 
along a channel in 3d, for which the constraint that the 
conflnement moves with the particle is not needed. 

We note that in the general case of a non-spherical par- 
ticle, or of a particle with non-uniform surface properties 
(e.g., a spatially varying effective mobility b) the particle 
can also rotate and the BC at the particle surface, Eq. 

should include a term accounting for a rigid-body 
rotation with angular velocity Jl. However, this angu- 
lar velocity turns out to be identically zero in the case 
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of a spherical particle with an effective mobility b which 
is constant on its surfaceiSi'^; therefore we completely 
disregard it here. 



3.3. Phoretic velocity 

The velocity V(t; 77) of the particle is obtained by re- 
quiring that the hydrodynamic force 



F = yy UerdS, 

\r\=R 



(7) 



where e^. is the radial unit vector and dS the surface 
area element on the spherical surface |r| = i? (note that 
we replaced R+ by R, as discussed in the previous sub- 
section), exerted by the fluid on the composite domain 
particle plus surface film vanishes [see the vector iden- 
tities in Eq. (jB7p . (i) - (iii) below this equation, and 
(r X (V X u))fj = xpdaUp — xpdpUa, with summation 
over P]: 



dS 



-p e,. 



+ ^V(ru) 

or r I r 



= 0. (8) 



\r\=R. 



If, as discussed above, a rotational motion with angular 
velocity ^l{t; rj) would have to be considered, too, this 
will be determined by the additional requirement that 
the motion is not only force free but also torque free^. 
This is again due to the fact that there are no net forces 
acting on the object composed of the particle and its 
surface film. 

The above argument for determining the velocity V 
has been discussed in detail by Anderson (sec Ref. [l^ 
and references therein), but it is often overlooked and re- 
placed by the incorrect argument of a balance between 
a drag force - i.e., the integral of the non-uniform os- 
motic pressure proportional to the density of solute (see 
Appendix [XJ - exerted on the particle and a Stokes- like 
viscous friction from the solvent (see, e.g., Refs.[3 l2^l23h . 
It is important to realize that the occurrence of composi- 
tion gradients in the solution does not give rise, by itself, 
to an osmotic pressure (see also Ref. [l5l ). in contrast to 
such an assumption made in Ref. [2^ . Such gradients will 
simply lead to diffusion of the product molecules, while 
the pressure in the solvent will adjust to accommodate 
the spatially varying chemical potential corresponding to 
the quasi-stationary density profile (mildly time depen- 
dent due to the overall increase of the solute number den- 
sity, in the case of the confined system), reflecting me- 
chanical equilibriuroi^. As a matter of fact, the origin of 
this osmotic pressure resides in the interaction between 
the particle and the product molecules, i.e., it requires 
the explicit consideration of the effective interaction po- 
tential between the particle and the product molecule 
(for a detailed illuminating discussion of this point see 
Refs. [T2llT7h . An intuitive argument that the use of such 



a "Stokes-formula" . which stems from a standard "drag 
balanced by viscous friction" type of reasoning, is incor- 
rect can be formulated as follows. At distances far from 
the particle the flow field should look like that produced 
by the superposition of a point force f at the origin, i.e., 
the center of the particle, which is the integrated (over 
the particle surface) effective "product molecules on par- 
ticle" interaction (the forces fp in Fig. [3] in Appendix 
E]) . and a distribution of effective point forces oriented 
radially (of the particle acting on the product molecules) 
in a small shell region around the surface of the particle 
(which is the aforementioned surface film), which upon 
integration gives exactly — f. Note that for repulsive effec- 
tive interactions and for an axisymmetric distribution of 
product molecules with an increasing density towards the 
source at z = —R, as in Fig. [H f = J^fxidS is oriented 
into the positive z-direction; for attractive interactions, f 
would be oriented into the negative z-direction. This can 
be seen as follows. For repulsive interactions, each of the 
forces f-p is oriented radially inwards; since the density 
of the product molecules is increasing towards the source 
located on the lower hemisphere, the magnitude and the 
projection onto the z-direction of the force due to the 
domain V located at any < 6* < 7r/2 is smaller than 
the one in the corresponding domain V located at n — 6. 
Thus the contribution into the positive z-direction from 
the lower hemisphere will dominate. For attractive inter- 
actions, the argument is simply reversed. This is in ac- 
cordance with the following two statements regarding the 
characteristics of the far-field (i.e., on length scales over 
which the particle plus surface film are seen as point-like) 
hydrodynamic flow: (i) In the absence of external body 
forces such as gravity or centrifugal forces the motion of 
the particle plus its surface film is net force free, (ii) 
The generated flow corresponds, within a first approxi- 
mation, to that produced by a "force dipole" [which is 
the superposition of a distribution of "force dipoles" (fp 
acting on the center and — f-p on the product molecules 
in 7)) as in Fig. [3] in Appendix [X] . or a higher order 
"force multipole", e.g., quadrupole (if it happens that 
the net force dipole vanishes, too)^ at the origin, i.e., 
the position of the center of the particle, rather than to 
the one due to a point force, which would be the case for 
an object uniformly dragged against the viscous Stokes 
friction (see, e.g., Ref. [241 ). For an unbounded system, 
these forces translate into a radial decay of the phoretic 
flow field proportional to (force dipole) or (force 
quadrupole), in contrast to the decay proportional to 
corresponding to a point force; for a bounded system, the 
differences between the flow fields cannot be any longer 
summarized by such a simple criterion as different power 
laws for the radial decay, but they remain significant nev- 
ertheless. These features of the hydrodynamic flow are 
discussed in more detail in the Appendix IB] (see also, c.f.. 
Fig. ED. 

Using the expansion of the velocity and pressure fields 
in terms of the solid harmonics K( (see Appendix [B|) . the 
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hydrodynamic force on the particle, defined by Eq. ([7]) 
and expressed as on the left-hand side of Eq. ([5]) , reduces 
to 



F = 47rp_2V[rPi(cos6')] , 



(9) 



where p-2 is the coefficient of A'_2 in the expansion of the 
pressurei^ [Eqs. (jB2|) and (jB3|l ]. and Pi is the Legendre 
polynomial of order one. All other terms vanish since the 
corresponding integrals are exactly equal to zero. Thus 
the requirement of a vanishing F implies p-2 = 0, which 
leads to (see Appendix [B)) 

IT 

b f 

V{t) =xi{v)- / de sin 6* cos 61 p{R,e,t;7]) (10) 



where 



XiW = 1 - 7^— 



2 7^5 



1 



(11) 



The structure of the expression on the right-hand side 
(rhs) of Eq. (fTUl) clarifies the meaning of the factor 
Xi(77), which varies between zero, at 77 — > 1, and one, 
at ?/ ^ 00 (see Fig. [2]). Without this factor, one has 
on the rhs the phoretic velocity in the unbounded space 
due to a source which generates a composition profile 
Poo such that PooJR, 0, t) = p{R, 0, t; rf) at all times t, as 
derived in Ref. [a [with A ^ — entering into 5; see 
the text around Eq. ([4])]. Thus xi(^) is a "hydrody- 
namic wall-correction" factor which quantifies and sum- 
marizes the effects solely due to the confinement induced 
changes in the characteristics of the solvent fiow. Note 
that Xiiv) < 1 means that the hydrodynamic effects due 
to confinement will tend to decrease the velocity from the 
value corresponding to an unbounded system. 



Density profile of tlie product molecules and 
diffusiophoretic velocity 



3.4. 



According to Eq. (jTU]), knowledge of the density 
p(|r| = R+,t;r]) of product molecules completely de- 
termines the velocity V of the particle as a function of 
r] < 00 for confined systems and thus allows one to quan- 
tify the effects of confinement on the resulting phoretic 
motion. 

Within the assumptions that the diffusion of product 
molecules is fast compared with the convection by the sol- 
vent flow, i.e., in the limit of small Peclet numbers, and 
that the product distribution p{r,t) is undisturbed by 
the flow, i.e., neglecting any so-called polarization effects 
of the surface filmi^, the time evolution of the number 
density p(r, t) of product molecules around the moving 
particle is governed, in the co-moving frame, by the dif- 
fusion equation 

dtp = DV^p + Bit)6{r-rs), R+ < \r\ < . (12) 



This equation is to be solved subject to the initial condi- 
tion (IC) of zero density of product molecules and to the 
boundary conditions of zero normal current at the con- 
fining wall (assuming that the latter is co-moved, with- 
out perturbing the solution, such that at all times it re- 
mains concentric with the particle) and at the outer edge 
r = of the surface film. Actually, the latter condi- 
tion should be imposed at the surface of the particle, but 
we shall consider only the case of a thin surface film so 
that due to its small volume the transport of solute to 
the surface film is negligible compared to the transport 
in the outer region and the zero normal current condition 
at R+ is a good approximation. Hence, the IC and BCs 



p(r,0) = 0. 



( dpir,t) 
\ dr 



0. (13) 



\r\=R+Mi 



We remark that the linearity of Eqs. (fT2|) and (fTO|) with 
respect to p imply that by superposition the solution of 
the present problem can be extended to the case of an 
arbitrary spatial arrangement of several catalytic sites on 
the particle surface. This will allow one to find an op- 
timal decoration for providing stability against particle 
rotations (caused by thermal noise) , which can otherwise 
spoil the unidirectional motion (see Ref. for such ex- 
amples of designed optimal surface distributions). 



Equation p2)) subject to the IC and BC conditions 
given by Eq. (|13p is solved by using the Laplace trans- 
form. Similarly to the approximation employed in Sec. 
13.21 from this point further in the calculations we replace 
everywhere i?+ by R because the difference S (the thick- 
ness of the surface film) between the two is negligible on 
the macroscopic length scales characterizing the trans- 
port in the outer region. The solution is obtained as 
a series in products of Legendre polynomials and modi- 
fied spherical Bessel functions of the first and third kind, 
respectively; the details are provided in Appendix [C] Fo- 
cusing on the particular case in which after its start the 
activity of the catalytic site is time independent, i.e., B(t) 
is constant: B{t) = H{t)/Tf, where H{t) is the Heaviside 
step function and Tf is the average production time for 
the creation of a product molecule. The Laplace trans- 
form of the integral in Eq. (jlOp is computed by using 
the Laplace transformed density p{r, C; 77) and the fact 
that cos6' = Pi (cos 6'). Formally inverting the Laplace 
transform, one obtains the diffusiophoretic velocity as a 
function of time and confinement: 



bxiiv) 



i?2' 7 TT^R-^DTf 



{c-'[mv)])\s, (14) 



where 



«l*3/2(yi) + (/3l-l)fc3/2(yC)l , (15) 
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^t:/{2z)Ii+iI2{z) and ^^+1/2(2:) 



\/ / {'2 z)K I j^i /2{z) are modified spherical Bessel func- 
tions of the first and third kind, respectively^^, and the 
dimensionless coefficients q;i(v^, '?) and /3i(>/Ci'7) are 
fixed by imposing the boundary conditions [see Eq. (jC8|) 
and Appendix [C] with ^ = C,R^ / D so that ^, and 
arc dimensionless]. 



4. DISCUSSION 

The complicated structure of the function 77) 
makes it rather laborious to carry out the inverse Laplace 
transform, so that the full time dependence of the veloc- 
ity V{t) cannot be derived easily. However, the asymp- 
totic (s ^ tD/R"^ > 1) value of the velocity V'^°°^r]) := 
lim V{s; rj) can be straightforwardly derived by using the 

s — >oo 

inversion formula for the Laplace transform^S, by notic- 
ing that has a simple pole at ^ = 0, which deter- 
mines the asymptotic value of the velocity as the residue 
of at ^ = 0^. The existence of such a constant 
asymptotic value, which at first glance seems to be in con- 
fiict with the fact that the density p(r; t) does not reach 
a steady state (since we consider a closed system with a 
source continuously producing non-interacting point-like 
particles), is due to the fact that the phorctic motion is 
determined solely by the gradient of the number density 
of the product molecules along the surface of the parti- 
cle. At long times, the total density of product molecules 
is large and any redistribution (which would lead to a 
change of the density gradient) proceeds on very slow 
time scales. (Ultimately the density of product molecules 
becomes so high that they no longer behave as an ideal 
gas and the latter assumption breaks down.) This yields 



where Vq 



n^R^DTf 
b 



Res^=o[mv)]=Vox{v), (16) 
y(°°)(7y 00) is the asymp- 



AirK^DTf 

totic velocity in the case of an unbounded system and 
[seeEq. (HH] 



Xiiv)X2iv), yv>l- (17) 



Thus X2iv) = {if + 2)/(?/^ - 1) is a "diffusion wall- 
correction" factor, and xiv) quantifies the combined wall 
effects, which in the present case factorize into a hy- 
drodynamic and a diffusion contribution. These two 
contributions oppose each other such that the confine- 
ment in hydrodynamics decreases the particle velocity 
(Xi < 1) whereas the confinement of the diffusion en- 
hances it (x2 > 1) (see Fig. [2]); the latter dominates so 
that X > 1 a-nd there is an overall enhancement of the 
velocity (see Fig. [5]). 

Assuming the Stokcs-Einstein relation for the diffusion 




FIG. 2: The long-time asymptotic velocity V^°°\'n)/Vo = 
xiv) > 1 (solid line) as a function of the confinement ry (see 
Fig. [l| and the liydrodynamic and diffusion "wall-correction" 
factors X'i-iv) < 1 (dotted line) and x^iv) > 1 (dashed line). 



coefficient of the product molecules in the solvent (see 
Sec. [1]) and using the expression Eq. ([3]) for the mobility 
6, the velocity Vq [Eq. can be rewritten as 



Vo 



3 A_ 



a 



(18) 



Surprisingly, this expression has no explicit dependence 
on the viscosity of the solvent or on the temperature. In 
fact, the dependence on viscosity drops out upon using 
the Stokcs-Einstein relation. In turn, the dependence on 
T is hidden in the temperature dependence of the (ef- 
fective) interaction between the particle and the prod- 
uct molecules which determines the thickness of the sur- 
face film, i.e., the length scale A [see Eq. (|4])]. In the 
case of only steric repulsion between the particle and the 
product molecules, for a particle of radius R ~ I fim 
[lOOnm], a product molecule with diameter a = 1.0 nm, 
A = 1.0 nm and a reaction rate (following Ref. @) of 
1/Tf = 25 kHz, one finds the velocity Vq — 10 pm/s [1 
nm/s] ; consequently V^°°^ (77) is of the same order of mag- 
nitude. Note that these values are much smaller than 
the estimate Vq = 0{iim/s) of Ref. H, which is based 
on a significantly larger value A^ ~ 10~^^ m^ than the 
one A^ ~ 10~^^ m^ corresponding to the reasonable esti- 
mate a = I nm for the diameter of the product molecule. 
However, this estimate corresponds to the case of a sin- 
gle catalytic site. One can argue that several such single 
reaction sites distributed closely around may lead to a 
significant increase of the velocity (intuitively, by a fac- 
tor equal to the surface density of such reaction sites). 
Thus the resulting velocities may eventually be closer to 
the experimentally reported values^ii^ Vq = 0(^m/s) for 
objects with a fim^ area covered by catalyst if the surface 
density of reaction sites is of the order of nm~^. More- 
over, if the effective interaction between the particle and 
the product molecules is attractive, the ratio (A/a)^ can 
be very largoi^. Another possibility leading to a signif- 
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icant increase is that of a classical hydrodynamic slip 
boundary condition with a large slip-length at the sur- 
face of the particle for the flow within the surface film, 
assuming that (see Appendix[A| classical hydrodynamics 
provides an accurate description for the solvent flow even 
at such small scales^^. (According to the text following 
Eq. (|A5|) . for the derivation of the effective mobility b 
in the phorctic slip- velocity [Eqs. ([1]), and (jA6p ] we 
have used a no-slip boundary condition for the flow in 
the surface film at the particle surface. This condition 
can be generalized to a slip condition characterized by a 
slip length.) Such hydrodynamic slip may be, e.g., due 
to roughness of the particle surface, and its efi^ect on the 
phoretic motion can be intuitively understood as follows. 
The slip over the surface of the particle facilitates the sol- 
vent flow within the surface film as compared to the case 
of zero slip (a "lubrication" effect), and thus for the same 
pressure gradient along the surface as in the case of zero 
slip the flow velocity in the surface film will be increased 
for nonzero slip- lengths. The phoretic slip Vg, which is 
the flow velocity at the outer boundary of the surface 
film, is thus also increased compared to its value in the 
case studied here so far that there is a no-slip condition 
at the surface of the particle; this leads to an increase 
(by up to orders of magnitude^^) of the phoretic velocity 
of the particle. 

The velocity V''°°\ri) of the particle can be positive 
or negative, depending on the details of the interaction 
potential. This means that the particle can travel ei- 
ther following the gradient of solute density, or against 
it. According to Eqs. pB]) and as well as the defl- 
nition of the effective mobility b, the sign of the phoretic 
velocity is opposite to the sign of the parameter b (or, 
equivalently, that of the parameter A), which in turn de- 
pends, inter alia, on the attractive or repulsive character 
of the effective interactions between the particle and the 
solute molecules [see Eq. Qj]. Equations (P) and ^ 
show that the sign of the phoretic slip velocity, i.e., of 
the polar component • eg is also opposite to the one 
of 6, because for our system dep is positive (sec Fig. [T|). 
In the case of hard core interactions only one has A < 
and thus 6 < 0; this implies V^°°^r]) > 0, i.e., the parti- 
cle moves in the direction of and thus away from the 
catalytic site, while the slip velocity points in the same 
direction as eg, which has a negative z-component, so 
that the flow around the particle is also oriented towards 
the catalytic sitei^. The confinement does not change 
the direction of V, which for an unbounded system is, as 
expected, in agreement with Ref. [^. A simple intuitive 
explanation for the direction of the motion of the parti- 
cle follows from the observation that there is no net force 
acting on the composite consisting of the particle plus the 
surface film, so that there is conservation of momentum. 
Because the solvent fiow around the particle is in the 
negative z-direction [see Fig. SUa)], the particle should 
move in the positive z-direction. Note that while this 
simple argument is clear in the unbounded case, for the 
confined system it breaks down because the formation of 



vortices [see Fig. SJ^c)] implies that there are spatial re- 
gions where the flow is in the positive z-direction. Thus 
a quantitative analysis is required. 

If 6 > 0, on the contrary, y(°°)(?7) < 0, which means 
that the particle will move towards the source of the prod- 
uct molecules. This is in contrast to the results presented 
in Ref. [2^ , which predict motion always away from the 
source. As discussed before the reason for this discrep- 
ancy is the ignorance in Ref. [13 of the mediating role 
of the solvent and its assumption of a mapping of the 
non-uniform density of the solute onto a non- uniform "os- 
motic pressure" acting on the particle (i.e., assuming that 
more product molecules impinge on the particle from the 
side with higher density), rather than mapping onto gra- 
dients in the solution pressure along the surface of the 
particle within its surface film. This opposing direction- 
ality in the case 6 > thus provides a simple criterion for 
an experimental discrimination between the predictions 
of the two proposals. 

Spatial confinement leads to an overall enhancement 
of the phoretic motion, because |Vb| < fo^' 
finite rj. As shown in Fig.[2l even at moderate values of rj 
the effect of the confinement is important, e.g., at 77 ~ 2 
there is a ca. 8% increase in the velocity compared with 
the unbounded case, and at ~ 1.5 the increase reaches 
ca. 25%. For 77 ^ 1 the velocity stays finite. This is, of 
course, an artifact which stems from the assumption that 
the product molecules are point-like. Actually, there is 
a lower cut-off i?f , where — R is of the order of the 
hard core diameter a of the product molecules, below 
which this assumption breaks down and Eq. (|16p is no 
longer valid. However, for a <C i?, which is a reasonable 
assumption, one has rjc ~ 1 + a/R > 1, and thus the 
steep increase of the velocity near = 1 is physically rel- 
evant. This is a clear counterexample for the statement 
in the conclusions of Ref. [13 that "the mobility of such 
swimmers will be hindered by channel boundaries" . 

Note that the velocity remains non-zero and finite for 
all values of ry, thus the two opposing effects described 
by the factors xi ^-nd X2 arc of similar magnitude. This 
can be understood qualitatively from the behavior of the 
velocity and density fields in the unbounded system: the 
hydrodynamic flow velocity decays as where r is the 
distance from the center of the particle^, and the gra- 
dients in the number density of product particles (for 
phoresis the gradients are relevant, not the density in 
itself) are are also varying as^ r^^ . Accordingly the 
confinement becomes relevant for both fields at simi- 
lar length scales and with a similar power-law behavior. 
Therefore one can expect that the confinement has an 
equally strong influence on the hydrodynamics and the 
diffusion. 



5. SUMMARY AND CONCLUSIONS 

In summary, we have discussed the effect on the 
phorctic velocity of a self-propelled particle due to a con- 
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fining wall for the solvent and for the reaction product 
emerging from a catalytic site on the surface of the dis- 
solved particle. 

The analysis is based on considering the present model 
in the context of the standard theory of phoresis. We 
have critically analyzed the assumptions involved in such 
an approach. Some of them are already present in the 
classical theory, others arise as a result of mapping the 
description of such "active" surface particles onto the 
framework of a theory developed to describe the case 
of inert particles immersed in a pre-defined, externally 
controlled concentration gradient. For example, within 
the standard theory of phoresis the perturbation of the 
steady-state, externally controlled concentration gradi- 
ent is induced by the interaction with the surface of the 
particle upon immersion. While this is an acceptable con- 
cept in that context, it is not clear how one can justify 
this if the concentration gradient develops as a function 
of time with the particle already present, as it is the case 
for a self-propelled object as the one we have discussed. It 
is beyond the scope of the present work trying to improve 
the general theory, but we consider a clear understand- 
ing of its limitations and of its possible shortcomings as 
a crucial step for both avoiding confusions such as those 
involving the application of a Stokes-forcc argument and 
for future theoretical developments. The main conclu- 
sion from this part of our work is that the development 
of a microscopic model for the dynamics in the surface 
film, eventually along the lines of Ref. which treats the 
solvent and the solute molecules on equal footing, seems 
to be very important. 

Within the confines of the standard theory of phore- 
sis, we have shown that the presence of a confining wall 
for the solvent and for the reaction product emerging 
from a catalytic site on the surface of a dissolved parti- 
cle leads to a significant increase of the velocity of the 
self-propelled particle. This results from two competing 
effects: an increase of the solute density gradients along 
the surface of the particle and a simultaneous increase 
of the hydrodynamic viscous friction. The former one 
dominates. If only steric repulsion between the particle 
and the product molecules is present, the absolute value 
of the velocity is expected to remain, in general, rather 
small. A direct experimental realization of the co-moving 
geometry considered here seems to be difficult. But the 
results which we have derived are expected to be appli- 
cable (at least qualitatively) for more general geometries 
(see Ref. [20| ). such as the motion of spherical particles 
along cylindrical tubes. In this sense, an experimental 
test of our results may be possible. 

Further extensions may focus on new phenomena 
emerging from more complicated geometries, as well as 
from relaxing the assumption of no interaction between 
the product molecules. By taking into account the actual 
finite size of the product molecules and by considering 
generic chemical reactions, in which the reaction prod- 
ucts are indeed different from the solvent molecules, it 
is expected that the production of such solute particles 



is accompanied by a non-uniform depletion of the sol- 
vent around the particle. Thus in general the motion of 
the particle will be determined by the gradients in both 
the solvent and the solute (see also Ref. 0). Moreover, 
such a depletion zone would also lead to a decrease of the 
production rate of the catalytic site and consequently, to 
a reduction of the density gradients. Accordingly one 
might expect that there are optimal values for the reac- 
tion rates for which the phoretic velocity is maximal. 

The interaction between the product molecules may 
play a significant role, and this deserves further discus- 
sion. If the density of the product particles is not low 
(which is reasonable to expect at least near the catalytic 
reaction sites, for fast reactions, and for slow diffusion 
of product particles) and the (solvent-mediated) solute- 
solute interactions have attractive or repulsive compo- 
nents longer ranged than the hard-core interaction dis- 
cussed above, several other effects have to be carefully 
considered, especially if they are as important as the 
particle-solute interactions. In this case, the distribu- 
tion of the solute particles in the direction normal to 
the surface of the big particle is no longer given by Eq. 
(|A3|) , but it is determined both by the particle-solute in- 
teractions and by the contributions of the solute-solute 
interactions, the latter depending on the whole distribu- 
tion of solute particles around the big particle. Thus it 
becomes a complicated non-local problem which cannot 
be easily addressed analytically. Within a mean-field ap- 
proximation, one may still assume that Eq. (|A3p holds if 
the potential is modified to include an averaged, effective 
solute-solute interaction. It is evident that even such a 
simplistic approach will lead to a different expression for 
the phoretic-slip velocity, and thus one can reasonably ex- 
pect qualitative differences compared to the predictions 
of the classic theory of phoresis. Moreover, the solute- 
solute interactions will lead to a dynamics which differs 
from simple diffusion and is determined, roughly speak- 
ing, by the density- and interaction-dependent collec- 
tive diffusion coefficient. For example, attractive solute- 
solute interactions will lead to a tendency of "clustering" 
and thus will hinder the relaxation of the solute gradi- 
ents, which intuitively would lead to an increase in the 
velocity of the particle, while repulsive interactions will 
help to dissipate the particle gradients and, intuitively, 
would lead to a decrease in the velocity. However, such 
intuitive arguments for changes in the velocity have to be 
carefully considered because, e.g., it is not clear if in the 
presence of significant solute-solute interactions the as- 
sumption of an ultrathin "unpolarized" surface film still 
holds. 
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APPENDIX A: CALCULATION OF THE 
PHORETIC SLIP VELOCITY 

The following considerations are connected to the cor- 
responding ones in Rcf. [l^ and put them into the context 
of the microscopic model discussed in Sect. [51 

Due to the azimuthal symmetry of the system, the flow 
field u(r) of the solvent has non-zero components only 
along the radial and polar directions. According to the 
considerations in Sect. [H the dissolved product molecules 
(solute) of diameter a are exposed to an effective inter- 
action with the particle within a surface film of thickness 
(5 ~ a. Accordingly, the solute molecules do not interact 
with the particle beyond the radial distance R+ = R + 5 
(measured from the center O of the particle). Since the 
particle radius R is much larger than a so that S/R <^ 1, 
this surface film of thickness 6 can be approximated to be 
locally planar (see Fig. [3]). In a small domain V centered 
at r = {R,9,(j)) (Fig. [S]), the local coordinate system (") 
is chosen such that x is along the polar direction eg while 
y is along the radial (normal) direction e^. 

Within the planar film approximation, which is reliable 
because the film thickness S is taken to be much smaller 
than the particle radius R {6 R), and for slow and 
smooth variations of the solute density along the surface 
of the particle (i.e., the solute density is assumed to vary 
over length scales of the order of the particle radius i?), 
the flow of the solvent relative to the particle surface is 
u{x,y) ~ {uj;{x,y),0) (with the meaning that the ratio 
Uy/iix is of the order 0{5/R) <^ 1)^. This can be in- 
tuitively understood starting from the incompressibility 
condition Vu = 0, where V indicates that the derivatives 
are taken with respect to the local coordinates. The com- 
ponent Ux is of the order of the slip velocity |vs|, which 
is expected to be proportional to the gradient along the 
surface, i.e., along the x direction, of the solute density 
p evaluated at Noting that dx = Rd6, one obtains 
dxUx ^ {1/R^)dgp. On the other hand, the component 
iiy has to rise steeply within the surface fllm of thickness 
6 from the value zero at the surface of the particle to 
the value corresponding to the radial component of the 
outer flow, which is expected to be also of the order of 
|vs|. Therefore dyiiy ~ |vs|/(5 ^ [l/{RS)]dgp. Since the 
derivatives of the solute density profiles with respect to 9 
are expected to not depend on the surface film thickness S 
and to be finite (due to assumed smooth, slow variations 
over length scales of the order of the particle radius), it 
follows that dyUy / dxUx ^ R/5 so that |9j}Uj}| ^ j^au^j. 
Thus the incompressibility condition dxUx+dyiiy = im- 
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FIG. 3: Schematic expanded view of a domain T) of the con- 
tact region between the particle and the solution, far from the 
reaction site, approximated locally by a planar geometry. The 
effective interaction between the reaction products, shown as 
horizontally hatched circles, and the particle is described by 
an effective potential ^'(y); y is the normal distance from the 
surface of the particle, also schematically depicted. The ar- 
rows at the product molecules indicate the forces — V<E'(i/) 
acting on them, while the arrow at the particle center shows 
the force f-D on it, due to the same interaction with the prod- 
uct molecules, which is equal in magnitude and opposite to 
the sum of — V^(y) over the product molecules positions in 
D. Here the case is shown that the product molecules in 
exert a net attractive force on the particle. The "outer edge" 
R+ of the surface film beyond which the effective interactions 
between the particle and the product molecules are negligible 
is indicated by the horizontal dashed line. The solvent, within 
which the product molecules move, is taken to be a homoge- 
neous background (not shown). Actually, S — 7?+ — R (of the 
order of the product molecules diameter a) is considered to 
be smaller than indicated here. 



plies that the leading contribution ~ 1/ {R6) must vanish 
separately because it cannot be canceled by subleading 
contributions ~ 1/R^. Therefore one has dyiiy ~ (up 
to terms of relative order 0{S/R)). This implies that Uy 
does not depend on y; because Uy = on the surface 
of the particle (independently of x), due to the surface 
being impermeable, it follows that ~ (in the sense 
of correction terms ^ 5/R) everywhere in P. 

As indicated in Fig. [31 the effective interaction of the 
product molecules (solute) with the particle, character- 
ized by the effective interaction potential "^{y) (in the 
limit of low solute density many-body effects can be ig- 
nored so that 4* is independent of x), gives rise to a 
force — V^" acting on the solute molecules. Assuming 
that the solvent can be approximated as a continuum on 
the length scale of the solute, the no-slip condition at 
the surface of the solute molecules implies that a body 
force density — pV^* is transmitted to the solvent. (Note 
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that, as discussed also in the main text, the approxima- 
tion of the solvent as a continuum breaks down if the 
solute molecules have a size similar to that of the solvent 
molecules, as in the case of the experiments discussed in 
Refs. in such a case this transmission of the body 
force has to be considered as an assumption to be veri- 
fied a posteriori.) Evidently, such forces arc present only 
within the surface film, where the effective interaction po- 
tential is non-zero. Within these assumptions, the Stokes 
equations for the flow VL(x,y) ~ (uj(x,y),0) within the 
surface film take on the form: 



(Al) 



The equation corresponding to the x component of Eq. 
()Aip can be further simplified by noting (with an argu- 
ment similar to the one used in the paragraph above) 
that the partial derivative I^Im^I is much smaller than 
(by a factor of the order 0[{5/Rf]). Thus for the 
X component Eq. (|Aip leads to 



whereas the y component yields 



(A2a) 



(A2b) 



We further assume that within the very thin surface 
film the relaxation of the solute density profile along the 
direction normal to the surface of the self-propelled par- 
ticle towards a steady-state (zero net current, i.e., the 
diffusion current due to gradients in p along the y direc- 
tion is balanced by the convective current generated by 
the force-field — V^") is fast compared to the diffusional 
relaxation time of its gradient along the surfac o^^i^^ . In 
this case, within the surface film the number density of 
the solute, assumed to behave like an ideal gas, is given 
by the Boltzmann distribution corresponding to the ef- 
fective interaction ^'(y): 



pix,y)^pix,R+)e-^^<-y^ 



(A3) 



(The prefactor p{x,R+) reflects the fact that ^'(y S) 
becomes negligibly small; we emphasize that Eq. (|A3[) 
applies only within the surface film.) We note that this 
assumption can break down, especially in the immediate 
vicinity of the reaction site, if the reaction rate is high. 
In such a case one has to explicitly consider the coupled 
equations of mass transport for the solvent and solute 
within the surface film (see Ref. 0); as already men- 
tioned before, we consider here only situations in which 
Eq. dMI holds. 



Combining Eqs. (jA2bp and (|A3[) . integrating with re- 
spect to y, and using that the inner and outer solutions 
for the hydrodynamic flow should match smoothly at i?+ 
so that at R+ the pressure reaches its "outer" solution 
value Poutix,y = R+), one obtains the following expres- 



sion for the pressure field within the surface film: 

p{x,y) =Pout{x,R+) + kBT[p{x,y) - p{x,R+)] . (A4) 

This is the so-called "osmotic equilibrium" condition (be- 
cause on the right hand side the term after the plus sign 
has the form known as "osmotic pressure" ) along the sur- 
face normal, which varies along the surface of the particle 
due to the gradient in the solute densit y^^i^^ . Note that 
Eq. (|A4p also applies only within the surface film. 



By: (i) combining Eqs. (|A2a|) . ((M)) . and (jXi)) . 
(ii) noting that dxPout{x, R+) can be neglected because 
dxPout{x, R+) ~ ^|vs|/i?^ (this is because the outer solu- 
tion satisfies the force-free (^' = 0) version of Eq. (|A1[) : 
by evaluating it at R+ and along the surface of the par- 
ticle and by noting that the slip velocity varies over the 
macroscopic length scale R, the above conclusion follows) 
while with \uxiy = ^-i-)| = |vs| and Uxiy = 0) = one 
finds from Eq. (|A2a|) that dxp{x,R+) ^ p\vs\/S'^, (in) 
introducing 



h{x,y) 



UbT dp{x, i?_|_) 



dx 



(iv) integrating once, and (v) using dyUx\R^ — (which 
holds because the inhomogeneity causing the flow van- 
ishes at the outer edge of the surface film), one obtains 



dyUx 



du h(x, u) —: ~g{x, y) . 



(A5) 



Integrating by parts, using the boundary condition 
Ux{y = 0) = of no-slip at the particle surface, and 
noting that ug{x,u) — > for u — > (because h(x,u) is 
bounded with respect to u, so that g{x, u) is finite for all 
u) and g{x, = [by definition oi g{x, y)], one obtains 
the phoretic slip-velocity Vs = Ux{x ~ -/?+) eg as 



Vs(^) 



g / duu h{x ^ R^6,u) . 



(A6) 



Since h(x, u) decays rapidly for u > the integral in 
Eq. (|A6p can be extended to infinity without causing a 
significant error: 



v.(0) 



ksT dp{x, R^ 



dx 



du \ 



-/3*(ti) 



-1] 



-ee b 



dp{x = R+e,R+) 
dx 



(A7) 



The derivation of the expression for the phoretic slip- 

1 • • ■ dp 1 dp 

velocity IS concluded by notmg that — = — — — corre- 
cts /t+ du 

sponds to the gradient of the solute number density along 
the particle surface [see Eqs. ([T])-(|31)]. 
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APPENDIX B: CALCULATION OF THE 
HYDRODYNAMIC FLOW IN THE OUTER 
REGION (BEYOND THE SURFACE FILM) 
AND OF THE HYDRODYNAMIC FORCE 



A general solution for the three-dimensional steady- 
state force-free Stokes equations, 



Vu = 0, 



(Bla) 



(Bib) 



in spherical coordinates has been obtained by Lambi^ii^. 
Here we briefly present the derivation, adapted to the 
present case which has an additional azimuthal symme- 
try, and we also compute the hydrodynamic force in the 
case of a flow subject to the boundary conditions given 
by Eq. del). 

The general solution is written as u = u/iom + Up, 
where Uhom is the solution of the homogeneous equations, 
i.e., /iV^u/iom — (so that the components of Uhom arc 
harmonic functions, i.e., they obey the Laplace equation) 
and Vu/io)ri = 0, while Up is a particular solution of Eqs. 
(|BTa)) and ([Bib]) . 



General solution in spherical coordinates in the 
presence of azimuthal symmetry 



The calculation of u proceeds by separately determin- 
ing the components Uhom and Up defined above. The con- 
struction of the component Up (r, 9) [a particular solution 
of the inhomogeneous Eq. (jBlaP ] starts from the obser- 
vation that the pressure field is also a harmonic function, 
i.e., V^p = [this follows from taking the divergence of 
Eq. (jBlaP J. Therefore, it can be expanded in terms of 
the solid harmonicsi^, which are the eigenfunctions of the 
Laplace operator in 3d: 



where 



with 



Kilr,e) = r^Qt{cose) J £ Z, 



Pi{cos9), for ^ > 0, 
P|c|_i(cos6l), for e<0, 



(B2) 



(B3) 



(B4) 



and Pe{cos6) is the Legendre polynomial of degree £. 

In this representation, a particular solution Up is given 
byl^^ 



Up = ^ Pi[Aer'^WKe + Ber Ke] 



(B5) 



where 

A, 



' + 3 



2e 



(B6) 



This can be checked by inserting Eqs. (jB5j) and (|B6p into 
Eqs. (jBlb . b) and by noting that Vr = 3 and rdrKe = 
IKg. Note that the pole in at ^ = — 1 is cancelled in 
Eq. (IBSI) . 

The construction of the solution Uhom starts from the 
identity 

V(ru,,<„„) = r X (V X u/^om) + Uho,n X (V X r) 
+ (u/iom • V) r + (r • V) Uh orn ■ 

(B7) 

Expressing the vorticity V x u/iom as VT [this is possible 



because V x (V x Uho 



V(Vuw) - V^Uho 



0] 



and denoting the product r • u/iojn by <i>, and noting that: 
(i) V X r = 0; (ii) (u,,o,„- V)r = Uhom, (iii) (r- V)u,,om = 
rdrUhom; (iv) r x (VT) = -V x (rT), Eq. (|B7| can be 
re-written as 



rdrUh, 



V$ + V X (rT) , 



(B8) 



Due to V X Uhom = VT one has V^T = so that T is 
harmonic, and due to V^<1> = 2Vuftom + r(V^Ufto„j) ~ 
also $ is a harmonic function. By inserting the series ex- 
pansions in solid harmonics [compare Eq. (jB2[) ] of u;iom, 
T, and $ into Eq. (jBSp and by using rdrKg — £Kg, one 
finds that the solution Uhom is given by 



= ^[t,V X (rX,) + l>,Vif,] 



(B9) 



where T^ and are the corresponding expansion coeffi- 
cients of T and $ as in Eq. (|B2p . 

The sum of Eqs. (|B5P and (|B9p provides the general 
solution u = Up + Ufiom expressed as a series expansion in 
spherical harmonicsi^ii^ with the coefficients {p^, T^, 
determined by boundary conditions on spherical surfaces. 
Note that for a problem without azimuthal symmetry the 
derivation of the solution proceeds analogously and the 
only change is that the Legendre polynomials Pi(cos9) 
{£ > 0) are replaced everywhere by a linear combination 
(with respect to the index —£< m < £) of the spherical 
harmonics Yim{9, (j)). 



2. Boundary conditions for problems with 
spherical symmetry 

The general solution derived in the previous subsection 
allows one to determine the hydrodynamic flow obeying 
Eqs. (jBlb . b) in a spherical domain with prescribed ve- 
locity on the boundaries of the domain (or at infinity). 
This is carried out in the usual manner by expanding the 
prescribed surface fields in terms of Legendre polynomi- 
als and then equating them with the series representation 
of the general solution evaluated at the boundary to de- 
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tcrminc the expansion coefficients {pi, Ti, However, 
the prescribed boundary conditions for velocity fields can 
be used to obtain an equivalent set of boundary condi- 
tions, which exploits the simplicity of the derivatives of 
the solid harmonics with respect to the radial coordinate 
and thus significantly simplifies the algebra^^; here we 
follow this approach, and discuss the most general case, 
i.e., without azimuthal symmetry. 



Let U denote the prescribed velocity field on a spheri- 
cal boundary |r| = c, i.e., u(r ~ c, 6, cj)) ~ U(6', (f)); in our 
case U = V -I- Vs for c — R+ and U = for c = i?i . 
(i) By multiplying u(r — c, 6, (f>) = U(0, (j)) with e^, one 
obtains the radial component u,. of the flow field 



ir{r = c,0,<j>) ^ er -Vie,^) Ur{9,, 



(BIO) 



where Ur denotes the radial component of the velocity U 

prescribed at the spherical boundary. 

(ii) The radial component of the vorticity is given by 



• (V X u) 



1 



r sm ( 



(Bll) 



where ug and are the polar and azimuthal components 
of u, respectively. On the spherical boundary |r| = c, 
within the right hand side of Eq. (jBll[) and ug can 
be replaced by [/^ and Ug, respectively. Thus the radial 
component of the vorticity obeys 

[r • (V X u)]|,,=, = r- (V X U) 

= -^[dg{smeU^)-d^Ug]. (B12) 
sm 



(iii) Since Vu ~ 0, one has 

= rVu = rdrllr + '^Ur 

+ -^[dg{sin9ug) + dsu^] . (B13) 
smt^ 

Since drUr{9,(f>) = 0, ug{r — c,9,(j)) = Ug{0,(j)), and 
Uif,{r ~ c,9,(j)) = U^{9,cl)), Eq. (jB13p renders on |r| = c 
the boundary condition 



r drUr\ I r—c 



-r(VU) 



(B14) 



Determining Ur, r • (V x u), and rdrUr from the so- 
lution u = VLhom + Up [Eqs. (jBSp and (|B9p ]. evaluat- 
ing the expressions at r = c, and equating them to the 
rhs of Eqs. ((BTO)) . (|Bl2l) . and ((Bll)) . respectively, af- 
ter expanding the latter in terms of spherical harmon- 
ics, yields the equations of condition for the coefficients 
{pi, , $^}. Applying this procedure for all boundaries, 
combining all the resultant conditions, and solving for 
the coefficients {pi, T^, $f } yields the velocity u in terms 
oi\J{9,(j3) = u(r = c,9,cj)). 



3. The outer hydrodynamic flow 

We now apply the general results discussed in the pre- 
vious subsections to the particular confined system shown 
in Fig. [U which exhibits azimuthal symmetry and obeys 
the boundary conditions in Eq. ^ on the spherical sur- 
faces |r| = i?+ ~ R and i?i = 7]R, respectively. Using the 
expansion of u = Up -I- Uhom in terms of solid harmonics 
[Eqs. (jBSP and (|B9[) ] and exploiting the azimuthal sym- 
metry (i.e., the solution does not depend on 0), splitting 
this series into two, corresponding to £ > and i < 0, re- 
spectively, changing in the latter the index of summation 
according to £ '—^ —{n + 1) , n> 0, and then renaming n 
by £, one obtainsi^ 



E 



2^(2^ + 3) 



r'+'Pi + 



1 



Pi{cos9), (B15a) 



•(V X u), = ^^(£+1) 



e>o 



r drUr = 



Pt{cos9), 
(B15b) 



£>0 



2^i{2£ + 3) 



2^(2^- l)r 



, (^+l)(^ + 2) ~ 



r-f-l-2 



(^+1) 



(cos 61). 
(B15c) 



Equations (|B15ap . (|B15bp . and (|B15c|) relate to Eqs. 
(|BT0)) . ((BT21) . and ((BTI)) . respectively. The rhs of Eqs. 
(|B15al) and (jB15cp do not depend on the coefficients {Tf } 
because the contributions T^V x {vKi) in Eq. (jB9p have 
no radial component. Similarly, the rhs of Eq. (|B15bp 
does not depend on {pi} and {$^} because r- (V x Up) = 
[Eq. jHH] and V x {VKi) = [Eq. (|B9])], respectively. 

By noting that VV = 0, V x V = [because the 
particle velocity V ~ Vgz depends only on t and 77, see 
Eq. ([TU)) ]. and V x = [because the slip velocity 
is given by a gradient, see Eq. ([T|)], one obtains that 
on the spherical surface |r| — Rj^ ~ R one has [see Eq. 
(|BT0)) ] Ur = {'V + Vs),. = Vcos9 [because V = Ve^ and 
e^Vs = see E g. jMI], r-[V xU] = r-[Vx(V+Vs)] = 
[see Eq. ((BT2I) ]. and [see Eq. fBTi)) ] 

- rVU = -rV(V + v,) = -rVv, = -k{9, R) , (B16) 

where [due to Vs = — 6Vsp(i?, 9; t, 77)] 

b 



k{9,R) 



■ dg {sin9dgp) 



Rsm9 

J2ki{R)Peicos9), (B17) 



i>0 
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p denotes p(R^6;t,r]), and 

2e + i 



kt{R) 



d9 sine k{e,R)PeXcost 



— I dO smP 



2 R 



de dO 



(B18) 



(The last equality follows upon integrating by parts.) By 
evaluating the right hand sides of Eqs. (|B15ap . (|B15b[) . 
and (jBlScp at R and Ri and equating them on the left 
hand sides with V cos{9),Q,—k{9,R) and with 0, 0, 0, 
respectively, the coefficients , $f} are obtained for 

any ^ > as the solution of the system of equations given 
by 



IR'+^ 



2p{2e + 3) 
+iR'^^^^e 



Pe 



eje + 1) R^+^ , 

2^(2^ + 3) 
£{£~l)R^-^^i 



-ie+1) = VSis , (B19a) 



= 0, (B19b) 



2p{2i-l)R 
{£+!){£ + 2) ~ 



jP-{i+i) 



Jli+2 



(«+l) 



—ki{R) , 
(B19c) 



£R{ 



e+i 



e + 1 



2p{2e + 3y^ 2p{2l-l)R[ 



JP-{I>+1) 



R\Te 



+1 



l)Rt 



e+i 



R\ 



£{t + l) 



0, 



(B19d) 



(B19e) 



2p{2£ + 3) 



2p{2l - l)R\ 



P-{i+i) 



+l[^£-\)R^ <l>fH -^j^ '±'-(£+1) 



R\ 



= 0. 
(B19f) 



Note that Eqs. (|B19bp and (|B19e|) imply that = 
T_f = for all I. Therefore, for any given I one is left 
with a linear system of four equations for four unknowns: 
P£, which (in the generic case) ad- 

mits a unique solution. 

Taking ^ = 1 in Eq. (jB19p . one obtains a system of 



four linear equations for the unknowns pi, $i,p_2, ^'-2, 
which depend parametrically on V, i?, and p, [see Eq. 
(|B6P ]. The absence of a body force acting on the particle 
implies p_2 [see Eqs. ([7]) - ©]. This leads to: 



Y 



kiR) 



5 7/ 

2r/5 



By noting that 



ki{R) 



/• ,„ . dpdPiicost 

/ do sni — 

2R / 961 de 



(B20) 



(B21) 



36 
'2R 



desin'e^ 

oO 



7T 

= 3— J de sine COS ep{R, 



one obtains the expression Eq. pUj) in the main text. 
Because dgp is positive for the system depicted in Fig. [1] 
the second equality in Eq. (jB2ip implies that the sign 
of fci(i?), and therefore that of V [see Eq. (|B20| ]. is 
opposite to the one of the effective mobility b. 

Before concluding this appendix, we note that the gen- 
eral solution for the hydrodynamic flow [Eqs. jBSj, (|B9p . 
and (jB19p ] allows one to gain insight into the qualitative 
differences between the phoretic motion mechanism and 
the one of an "osmotic pressure" propeller—. In Figs, 
mja) and (c) we show the normahzed velocity fields u/|u| 
corresponding to diffusiophoretic motion due to steric re- 
pulsion only. In this case 6 < [see Eq. (U)] so that the 
particle velocity V = Ve-^ is oriented along the positive 
z-direction, i.e., away from the source. Figure [4ja) cor- 
responds to an unbounded system whereas Fig. \^c) de- 
picts the corresponding confined system with 77 = 5; both 
cases refer to the same (chosen) constant surface gradient 
dgp{R,6;rj) > of the solute number density. (Being a 
constant, this gradient scales out and is absorbed in the 
velocity scale. Therefore the dimensionless ratio u/|u| 
is independent of it.) Since h enters the problem only 
via the multipHcation of dep{R,e;t,'ri) [see Eq. (|B2ip ]. 
also the absolute value of b drops out from the normal- 
ized flow field. The viscosity p drops out, too. For the 
unbounded case, the solution [see also Eq. (35) in Ref. 
nil is determined from Eqs. (|B19ap . (|B19bp . and (|B19cp 
and the requirement that the flow field vanishes at infin- 
ity. This implies that in Eq. (jBlSp all coefficients mul- 
tiplying terms ^ with £ > are zero. In the confined 
case, the flow field is approximated by keeping terms up 
to ^ = 50 [Eqs. (|B19p ] in the series expansion of Eqs. 
(|B5P and (|B9P : this provides a reasonably good approx- 
imation, except near the boundary at i?i and for polar 
angles close to zero or tt where apparently there are de- 
viations from the no-slip boundary condition. However, 
one should keep in mind that in this normalized rcprcsen- 
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FIG. 4: Normalized velocity fields u/|u| in the laboratory 
frame corresponding to the diffusiophoretic motion due to 
steric repulsion only [h < 0) between the solute and the parti- 
cle for the case of (a) an unbounded system and (c) a confined 
system (r; = 5) . The panels (b) and (d) show the correspond- 
ing results in the case that the motion occurs as a result of 
a drag provided by the "osmotic pressure" integrated along 
the surface of the particle. The dashed circles in (a) and (b) 
indicate the position at which the confining wall is placed in 
(c) and (d). In all cases (a) - (d) the particle moves upwards, 
i.e., its velocity V is in the positive z-direction (as shown by 
the vector on the particle) . Here we consider only time scales 
on which the movement of the particle does not change the 
fiow fields. 



tation of the flow u/|u| small numerical deviations from 
u — at i?i arc sharply overemphasized. For compari- 
son, in Figs, m^b) and (d) we present the corresponding 
results for the case in which the motion of the parti- 
cle occurs as a result of a drag provided by an "osmotic 
pressure" integrated over the surface of the particle. In 
this case the particle velocity is oriented along the posi- 
tive z-direction, too, because of the assumption that the 
pressure is proportional to the solute number density and 
thus is larger at the lower hemisphere. In this case Eq. 
(|B19[) is solved under the assumption of a no-slip bound- 
ary condition at R, i.e., all the coefflcients ki are set to 
zero; this allows one to obtain a solution in closed form 
in both the unbounded and the confined cases. Note that 
in this case p_2 is non-zero. It is fixed by the given drag 
[integrated "osmotic pressure", Eq. ©], and thus Eq. 
(|B19[) . evaluated at £ = 1, leads to a relation between 
the velocity V and the drag force. (In the unbounded 
case, this relation is the well-known Stokes formula for 
the viscous friction on a sphere.) In Eq. (|B19a[) V sets 
the velocity scale, which drops out (as well as the drag 
force, because V is proportional to it) of the dimension- 
less ratio u/|u|. 

Qualitative differences between these two mechanisms 



are particularly clearly visible in the case of the un- 
bounded systems. For example, the flow along the z- 
axis is in the negative direction (i.e., opposite to the 
motion of the particle) in the phoretic case, while it 
is in the positive direction for the "osmotic" propeller; 
the flow lines for the phoretic motion are closing back 
and then "slide" along the surface of the particle. Since 
closed-form analytical expressions are available in these 
unbounded casesi^*^, the qualitative differences noticed 
above can be correlated with the different characteristics 
of the flows. In the case of phoretic motion the flow field 
decays as at large distances r from the particle, which 
is a "multipole" far-field corresponding to the action of a 
force-quadrupole disturbance, in contrast to the de- 
cay in the osmotic case, which is a "monopole" far-field 
corresponding to the action of a point-force ( "Stokeslet" ) 
disturbance. The faster decay of the velocity field in 
the first case means that a boundary at Ri will perturb 
the phoretic flow much less. This explains why for the 
phoretic motion the viscous friction increases with the 
conflnement weaker than what would be expected from 
naively using a Stokes- friction concept. 

The flow fields remain qualitatively different even in 
the presence of the confinement (in particular with re- 
spect to the opposite direction of fiow along the z-axis) , 
although the differences are not as pronounced as in the 
unbounded case. A further qualitative difference between 
the two confined flows is that in (d) a vortex structure 
is formed fully detached from the particle whereas in (c) 
the vortex involves flow along the surface of the parti- 
cle. Due to the presence of the boundary at Ri, all the 
terms ~ in the expansions of Eq. (jB15p are present 
in the solution and an analysis of the flow in terms of 
fundamental solutions, as for the unbounded case, is no 
longer possible. Note that since V x Uhom = VT and 
T = [see the text following Eqs. ([13 and (|BT9^ ] one 
has V X u = V X Up so that with Eqs. (|B5p and (|B6|) 
one obtains 



V X u = — e 



^ pe dK,[r,e) 



(B22) 



which is in general nonzero [e.g., for "osmotic" flowp„2 7^ 
(Eq. dS]))]. Thus in agreement with Fig. |3]the flow field 
can contain vortices. 



APPENDIX C: CALCULATION OF THE 
NUMBER DENSITY OF PRODUCT 
MOLECULES 

Equation (|12p . subject of the IC and BC conditions 
given in Eq. (|13p is solved by using the Laplace transform 
and the inversion theorem^^: 



/(C)^£[/]= dte-^'f{t) 



(CI) 
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and 



f{t)^C-'[m]:=— J dCe^'m, (C2) 
7—200 

respectively, where /(C) (the Laplace transformed quan- 
tities are indicated by an overbar) is assumed to be well- 
defined for ( S M+ and 7 g K is sufficiently large such 
that all singularities of /(C) lie to the left of the inte- 
gration path. By taking the Laplace transform of Eqs. 
(fT2|) and (fT3|) and by using the IC condition of zero num- 
ber density of the product molecules, one finds that the 
Laplace transformed density is given by 



p(r,C)=-a(C)G(r,g), 



(C3) 



where a{() = B{C,)/D and G{y, q) is the Green's function 
for the Helmholtz operator — q^, q ~ \JCJl5 > 0, 
satisfying the BCs of vanishing normal derivative at |r| — 
i?, i?i . 

Decomposing G{r,q) as G{r,q) = Gs(r, 9; r^) +.g(r, 9), 
where the singular part Gg is the free space Green's func- 
tion for the Helmholtz operator (which is known in any 
spatial dimension d, see, e.g., Ref. [28h . the initial prob- 
lem is reduced to that of finding the solution g of the 
homogeneous Helmholtz equation subject to the bound- 
ary conditions 



f dgjr, q) 
\ dr 



\r\=R.,Ri 



dGs{Y,q;Ys) 
dr 



(C4) 

In 3d the singular part of the Green's function is given 



q e 



An q\r - 



(C5) 



while the regular part (i.e., the general solution of the 
homogeneous Helmholtz equation) can be written as 

9{r,d,q) = ^[ae-ie+i/2{qr) + (3e ke+i/2{qr)]Pt{cos9), 

(C6) 

w here ie+ i/2{z) = ^/TT/(2z)Ie+i/2{z) and ^+1/2(2) = 
y^Tr/ {2z)Kg^i/2{z) are the modified spherical Bessel 
functions of the first and third kind, respectivelj*^^, 
is the Legendre polynomial of degree £, while the coef- 
ficients and (3^ will be fixed to fulfill the boundary 
conditions. 

Using one of the addition theorems for the Bessel 
functions2& and noting that the angle between r and 
is, IT — 9 (see Fig. [T]), Gg can be re- written as 



Gs{r,q;rg) = - ^ce{q)ii+i/2{qr<)ki+i/2{qry)Pe{cose) 

(C7) 



where Ci{q) = (-1)^(2^+ l)g/(27r2), r< = min(r,r,), 
and r> = max(r, Ts). By combining Eqs. (|C6i IC7I 



and the boundary conditions [Eq. (jC4[l ]. by noting that 
rs~R + e,e\,0 {i.e., the source of particles is on 
the surface of the particle), and by re- writing the coeffi- 
cients as ai = Ciii+i/2{qR)ae and Pt = ciii+i/2{qR)(3g, 
the dimensionless coefficients ai(qR) and f3g{qR) are de- 
termined by the solution of the following closed system 
of two linear equations with two unknowns (a^ and f3e): 

at + i3iWe+i/2iqR) = ve+i/2iqR), (C8a) 
ae + f3eWi+i/2{mR) = we+i/2{mR), (C8b) 

where ^£+1/2(2) := kf.+i/2{z)/ie+i/2{z) and We+i/2{z) := 
[dkg^i/2{z)/dz]/[dii^i/2{z)/dz]. The solution ai and 

Pi acquires a dependence on 77 via Eq. (jC8b[) : = 
ae{qR,T]) and Pe = Pe{qR,r]), with q_R = R^Q/D. 
This determines the Green's function G(r, q) and thus 
the Laplace transformed density p(r, C,) [Eq. (|C3p ] from 
which, in principle, the density p(r, t) is obtained via the 
inverse Laplace transformation. Note that although the 
final result is in the form of a series, the quantities in 
which we are interested, in particular the phoretic veloc- 
ity V{t) [Eq. [TU], involve integrals over the polar angle 9 
of the product between this series and a specific £ Legen- 
dre polynomial Pg (cos 9) [£ = \ in the case of the velocity, 
^ = in the case of the total number of product molecules 
(see below)] and thus only one of the terms from the se- 
ries will contribute. 

We make the following three remarks: 

(i) Since for qR > one has lim u'^+i/2('79^) — 0, in 

rj — 'oo 

the limit 77 ^ 00 one finds = and Eq. (|C8ap reduces 
to the corresponding BC in Ref. 0, i.e., as expected the 
solution for the unbounded case is recovered. 

(ii) Equations (|C8a[) and (|C8bp will not coincide in the 
limit 77 — > 1 because even in this limit of extreme confine- 
ment a source singularity is present at r,, and is picked 
upby G,(g) [Eq. jCB]. _ 

(iii) The Laplace transform N{p) of the total number of 
product particles in the system at time t is 

N{p^ Dq^) = 271 j drr^ j d9 sin9Po{cose)[-~a{p)]G{q) 

R. 

= -'iTTa{p)co{q)ii/2{qR) 
Rl 

X J drr^[aoti/2iqr) + {Po - l)h/2iqr)] (C9) 

R 

t 

= Da{p)/p =J> N{t) = J dt'A{t') , 


i.e., as expected N{t) is given by the time integral of the 
production rate, providing a welcome consistency check. 
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